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Abstract. This article is devoted to the analysis of eigenfunctions (modes) and approx- 
imate eigcnfunctions (quasi-modes) of the Laplacian on a compact manifold (M, g) with 
completely integrable geodesic flow. We give a new proof of the main result of [TZ] that 
(M,g) with integrable Laplacians and with uniformly bounded eigenfunctions must be flat. 
The proof is based on the use of Birkhoff normal forms and on a comparison of modes and 
quasimodes. In the process, we discuss tunnelling between resonant tori and give a proof 
that eigenfunctions concentrate on individual tori in the non-resonant case. We also give 
brief expositions of results in [TZ2, SZ|. 



0. Introduction 

In this article we are interested in the semiclassical h — > asymptotics of modes (eigen- 
functions) 4>j(h) 



(h 2 A + V% = Aj(/i)0 i5 (fa, cf> k ) = 5 jk (1) 

of Schroedinger operators on a Riemannian manifold (M,g). Usually we take (M,g) to be 
compact for convenience, but equally interesting problems occur on non-compact manifolds 
such as M n . Many questions about eigenfunctions are motivated by modern mesoscopic 
physics and quantum chaos, in particular by the physics of quantum dots, wires and corrals 
(for a recent sample, see e.g. WE , |LHrh\lLE| . [K[. [Mi| , Mi2 |). Let us mention just a few topics 
and recent articles from the mathematics and physics literature : 

• Norms of eigenfunctions: \\fa\\ L v ||ABST| , |A~S| [Bollr] , 0, || gZ[ ; 

• Eigenfunction scarring flAAj [Bell , 0, [RSI, [W]]; 

• Volume and distribution of nodal surfaces {fa = 0} PGS[ [BJJ |CMKNG| , pF| , pF2t |JL 



• Number and distribution of critical points {V<f)j = 0} ||JN| , |5Z| ; 

• Matrix elements (A(j)j,(j)k) of pseudodifferential operators |[KH| , |Z[]; 

• Resonance and tunnelling between wells or tori |[HS| , pSuj , |M] , |W| . 

Most of the classical and more recent problems about eigenfunctions are believed to be very 
difficult. The main reason for the difficulty is that in general there are no explicit formulae 
for eigenfunctions or even for approximations to them. Lacking such explicit formulae, one 
has to rely on local approximations that lose track of the global geometry of (M , g) or on 
wave equation or resolvent methods, which involve averages over many eigenfunctions. It is 
thus difficult to obtain global information about individual modes. 
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In this article, we would like to review and revisit some recent results about L p -norms 
of eigenfunctions, obtained in part by the authors and in part by one of the authors in 
collaboration with C. Sogge. For some (M,g), L p - norms of eigenfunctions are extremally 
large ||SZ|| and for some they are extremally small ||TZ| , [TZ2|| . The problems studied in 



||TZ| , [rZ2 , were to characterise metrics and potentials producing extremal growth rates 



of eigenfunctions. Of course, one would like to have some control over growth rates of 
eigenfunctions for any (M, g, V). Our main focus in this article is to give a different and 
in some ways more natural proof of the results on L°° norms of modes in the quantum 
integrable case of [ |1'Z|| , based on the construction of quasimodes. 

On a more expository level, we would also like to advertise quantum integrable systems as 
a simple setting for harmonic analysis. Quantum integrable Laplacians are generalizations 
of Laplacians on flat tori and of harmonic oscillators, and have most of the features which 
make these settings attractive for harmonic analysis. Yet they also bring in interesting 
new geometric and analytic phenomena. Most of the special eigenfunctions studied in the 
literature are quantum integrable, see e.g. [[B.K.S| , [Bourj , |CV1| , |D|, [IN], IKMSf . To make 



the ideas more accessible to harmonic analysts, we describe a number of two-dimensional 
examples in some detail in Section (||). 

To motivate the problem and result, let us consider the setting where we first learn about 
eigenfunctions, namely the flat torus ~R n /L where L is lattice such as Z n . The eigenfunctions 
<p\ are of course the exponentials e^ A,a ^ where X E L*, the dual lattice. The harmless looking 
facts 

A ( x ) =e ^\ |0 A (x)| = l, (Vx) 

about the most familiar eigenfunctions raise the questions: when can eigenfunctions be 
expressed in this kind of WKB form a(x)e lS ^ x ^ h or as an oscillatory integral? And how 
common is it that eigenfunctions have uniformly bounded sup norms? 

Heuristic answers to these and other questions about eigenfunctions are suggested by 
correspondence principles between excited states of quantum systems and their underlying 
classical dynamics. The quantum Hamiltonian is simply the Laplacian A, while its classical 
counterpart is its bicharacteristic flow, namely the geodesic flow G l : S*M — > S*M. The 
clearest statements concern the most extreme dynamical regime, integrable systems and 
chaotic systems, with mixed and somewhat fuzzy statements about intermediate systems. 
The following description summarizes the current conjectural picture of eigenfunctions; it is 
not intended to be exhaustive. Eigenfunctions are described according to the regime of the 
corresonding classical dynamics. 

• The Integrable regime: Approximate eigenfunctions may be constructed by a WKB 
method associated to invariant tori satisfying Bohr-Sommerfeld quantization conditions. 
Sequences of eigenfunctions localize in phase space on (unions of ) the invariant tori for 
the geodesic flow. They 'scar' (concentrate singularly) on the projections of those tori. 

• The intermediate KAM regime: WKB quasimodes may be constructed for a positive 
measure set of "Cantori" . There are regular and irregular modes corresponding to the 
Cantori and the regions of instability between them. 
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• The chaotic regime: eigenfunctions behave like random waves, e.g. like random spher- 
ical harmonics. They do not localize in phase space, but rather become uniformly 
distributed, as do their critical points and zero sets. Only low-lying modes should scar. 

Quasi-modes are essentially parametrices for modes, in some ways analogous to parametri- 
ces for operator kernels (e.g. the Hadamard parametrix for the heat or wave kernel). Unlike 
operator parametrices, however, they do not necessarily approximate individual modes. We 
return to this point below. 

Quasi-modes are usually thought of as oscillatory integrals associated to geodesic-flow in- 
variant Lagrangean (or isotropic) submanifolds of an energy level {H(x, £) = + V(x) = 
E} C T*M. There are axioms for quasimodes | fj V2| | which allow for more general construc- 
tions, in particular in the KAM case [|CV2| , [[], [FJ. Existence of quasimodes associated to 
a Lagrangean submanifold A involves a Bohr-Sommerfeld quantization condition on A and 
some stability properties. For instance, there are always quasi-modes associated to stable 
closed geodesies |CV2j , [Rj, but there are no known quasi-mode constructions associated to 



unstable closed geodesies. For this reason, quasi-mode constructions are most effective in 
integrable or KAM systems, and these regimes represent the best chance of understanding 
individual eigenfunctions. In the opposite chaotic regime, there are no known quasi-modes 
and heuristic reasons suggest there aren't any explicitly constructible approximate solutions. 
Indeed, the modes are believed to behave like random waves and to be rather featureless. 
Finally, in the no man's land between the nearly integrable and the quite chaotic, there are 
few if any heuristics or results of a global nature. The only known results are obtained by 
piecing together local information on eigenfunctions. 

Let us now state some recent rigorous results. First we need some notation: by \l/ m (M) we 
denote the space of mth order pseudodifferential operators over M. We say that an operator 
P G # X (M) is quantum integrable if there exist P 2 , . . . , P m G ^ l (M) such that [Pj, P k ) = 
and such that Pf + • • • + P^ is elliptic. We denote by p k = (Jp k the principal symbol of P k . 
Since crp^Pj] = {PuPj} (the Poisson bracket), it follows that the p/s generate a homogeneous 
Hamiltonian action $j of t G W 1 on T*M — with moment map 

V:T*M-0^R n , V=( Pl ,...,p n ). 

We denote the image V{T*M — 0) by B, by B reg (resp. B sing ) the regular values (resp. 
singular values) of the moment map. 

We are going to impose a 'finite complexity' assumption on the classical integrable system. 
Suffice it to say that it holds for all systems of interest in physics: for each b = (b^\ b^) G 
B, let m c i(b) denote the number of M ra -orbits of the joint flow $ 4 on the level set V~ l (b). 
Then 

Finite complexity condition : 3M : m d (b) < M (V6 G B). (2) 

When b G B reg , then V~ l {b) is the union of m c i(b) isolated Lagrangian tori. If b G B sing , 
then P _1 (6) consists of a finite number of connected components, each of which is a finite 
union of orbits. These orbits may be Lagrangian tori, singular compact tori (ie. compact 
tori of dimension < n) , or non-compact orbits consisting of cylinders or planes. 
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On the quantum level we will assume: 

Bounded eigenvalue multiplicity : 3M' : m(A) < M' (VA; m(A) = dimV\). 

(3) 

This reduces estimates of eigenfunctions to estimates of joint eigenfunctions of the commuting 
operators. When unbounded multiplicities occur, such as on M n /1/ 1 , then it is well-known 
that one may construct sequences of eigenfunctions with unbounded sup norms even on flat 
tori. 

Theorem 0.1. Suppose that A is a quantum completely integrable Laplacian on a compact 
Riemannian manifold (M,g), and suppose that the corresponding moment map satisfies (^). 
Then: if the IS.- eigenfunctions have uniformly bounded L°° norms, then (M,g) is flat. 

More generally, suppose that h 2 A + V is a quantum completely integrable Schroedinger 
operator, and that the corresponding moment map V is proper and satisfies ^). Assume 
there exists an energy level E such that the eigenfunctions have uniformly bounded L°° norms 
as h — > 0. Then: if E > max V, and (M, (E — V)g) is flat. If (a) (or (b)) holds for all 
energy levels E in an interval E% < E < E 2 , then (M,g) is flat and V is constant. 

To prove the result, we must relate the analysis of modes to the geometry of the Liouville 



foliation (by invariant tori) of T*M of the completely integrable system. In |TZ|| , this relation 
was brought out through trace formulae. In this article, we approach it through quasimodes. 
With some work, we believe the results should generalize to some degree to quantum KAM 
systems, i.e. to small perturbations of completely integrable systems. However, we leave 
that extension for the future. 

Roughly speaking, the idea is this: Quasi-modes are directly associated to the individual 
Lagrangian tori of the completely integrable system. Hence we can directly relate norms of 
quasimodes to geometry of the tori. To draw conclusions about norms of modes, we must 
then relate modes and quasimodes. Since the classic article "Modes and Quasi-modes" 
of V. I. Arnol'd (also based on examples due to A. I. Schnirelman), attention has been drawn 
to the danger of of drawing conclusions about modes from conclusions about quasimodes. 
As Arnold cautioned, modes are not always approximated by quasimodes, but rather may be 
linear combinations of quasimodes with close quasi-eigenvalues. If the linear combinations 
involve 'many' closely spaced quasi-eigenvalues, then the connection between modes and 
quasi-modes becomes blurred. However, this danger mainly pertains to scalar spectra, i.e. 
the spectrum of one operator; it does not occur in the approximation of joint eigenfunctions 
of a completely integrable set of operators by joint quasimodes, since the joint spectrum in 
M. n is much better separated than the scalar spectrum of the Hamiltonian alone. For the 
joint spectrum, a simple non-resonance condition is sufficient to ensure that modes are even 
well-approximated by individual quasi-modes. 

When modes of integrable systems can be connected with invariant tori for the geodesic 
flow, one can hope to relate growth rates of L p -norms of sequences of eigenfunctions to 
singularities of projections of the associated tori. For instance, as will be discussed further 
in Section (|5[), invariant eigenfunctions on surfaces of revolution correspond to to invariant 
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torus in S*S 2 consisting of meridian geodesies between the poles. This torus has a "blow- 
down" singularity under the natural projection n : S*S 2 — > S 2 , i.e. the fiber S^ P S 2 — > {NP} 
over the north pole is blown-down to a point under the projection. This kind of singularity 
causes the largest concentration of the mass of a mode at a point; hence the invariant 
eigenfunctions which correspond to the meridian torus are (at least in order of magnitude) 
extremals for the L°° and other L p norms for sufficiently large p. Other important tori in 
the LP theory are the singular tori where the torus drops dimension, e.g. the equators of 
a surface of revolution. Eigenfunctions concentrated on such singular tori are (in order of 
magnitude) extremals for low L p norms: although they do not blow up as much, they blow 
up on a set of larger measure and hence are favorable for LP norm blow-up. 

Let us now recall some other recent results; we refer to ||TZ2|, [SZ|| for further discussion. 



The next result also pertains to quantum completely integrable systems. We impose a 
stronger hypothesis on the classical system, namely Eliasson non- degeneracy (see flTZ2| for 



background). We then obtain quantitative results on the LP blow-up of eigenfunctions in 
non-flat cases, i.e. in cases where (by the previous theorem) some blow-up must occur. We 
emphasize that our notation for eigenvalues ([!]) is such that A<pk = \\4>k- 

Theorem 0.2. Suppose that {M,g) is a compact Riemannian manifold with completely in- 
tegrable geodesic flow satisfying Eliasson's non- degeneracy condition. Then (M, g) is non-flat 
if and only if for every e > 0, there exists a sequence of L 2 -normalized eigenfunctions (pk 
satisfying: 

\<Pk\\oo>C(e)\r € . 



4p 



\\Mp>C(e)X k 

Modulo the e's, these estimates are sharp. They are based on the existence of codimension 
one singular leaves of the Lagrangian foliation, and on the construction of precise quasimodes 
attached to non-degenerate singular leaves. In general, a completely integrable system does 
not have higher codimension singular leaves (e.g. one may take the product M = T n ~ 2 x S 2 
of a flat torus and a surface of revolution to obtain completely integrable (M, g) of any 
dimension whose only singular leaves are of codimension one). When higher codimension 
leaves occur, the blow-up rate increases. In the case where (M, g) has a one- dimensional 
singular leaf , then our methods show that there exists a sequence of eigenfunctions satisfying: 

n-l 

||0fe||oo > C(e)X k 4 . Such a one-dimensional leaf is simply a stable elliptic geodesic 7 which 
is fixed by the other Hamiltonian flows of the completely integrable system. This estimate 
is comparable to the recent estimate of Donnelly ]B| on the blow-up rate of a sequence of 
eigenfunctions on certain (M, g) with an isometric S l action. Note that Donnelly denotes 
the eigenvalue by A while we denote it by A 2 . Donnelly also assumes the existence of a 
non-degenerate singular orbit 7 of the joint action of the geodesic flow and S 1 action. Such 
a 7 is a a stable closed geodesic, such as the equator on a convex surface of revolution, 
which is invariant under the S 1 action. The corresponding sequence of eigenfunctions is 
actually a sequence of singular quasi-modes associated to 7. Thus, it is the same phase 
space phenomenon in either setting which is responsible for this common blow-up rate. 



6 



JOHN A. TOTH AND STEVE ZELDITCH 



We now turn to results by one of the authors with C. Sogge on the opposite extreme: which 
(M, g) have maximal eigenfunction growth? It follows from classical remainder estimates of 
Avakumovic and Hormander that L 2 -normalized eigenfunctions A0>, = X 2 4>\ have L°°-norms 
satisfying 

H^aIIoo < C\—, n = dimM. 
When a sequence of eigenfunctions of (M, g) actually achieves this bound, we say that (M, g) 
has maximal eigenfunction growth. The main result of ||SZ|| implies a necessary condition on 
a compact Riemannian manifolds (M, g) with maximal eigenfunction growth: there must 
exist a point x G M for which the set 

£ x = {teS:M:3T:e Wx TZ = x} (4) 

of directions of geodesic loops at x has positive surface measure. Here, exp is the exponential 
map, and the measure \Q\ of a set Q is the one induced by the metric g x on T*M. For instance, 
the poles Xn,%s of a surface of revolution {S 2 ,g) satisfy \C X \ = 2n. More generally, we may 
consider maximal IP growth of eigenfunctions. The universal esimates in this case are due 
to Sogge Pol , |So2|[ . 



Theorem 0.3. If \C X \ = for every x e M then 



^ = (A'«), P > 2 -^ (5) 

L 2 (M) 



where 



2(n+l) 



< p < OO 



^ 2 V2 p)l Z — V — n-1 ■ 

It is proved that generic metrics g on any M satisfy the hypothesis of having zero-measure 
loop directions at every point. Hence, maximal eigenfunction growth is a rare condition. 
The strongest conclusion is available for real analytic Riemannian manifolds: 

Theorem 0.4. Suppose that (M, g) is a real analytic with maximal eigenfunction growth. 
Then there exists a point m G M such that all geodesies issuing from the point m return to 
m at the same time t. In particular, if dim M = 2, then M is topologically a 2-sphere S 2 
or a real projective plane M.P 2 . 

It would be interesting to investigate maximal eigenfunction growth in small IP norms. 
It is explained in [5ol , |So2|| that the geometry of eigenfunctions changes at the critical 



V = 2 ^n-P ' an d si zes °f 1°°P directions is no longer relevant. In the case of the found 
sphere, at this value of p the extremals for the LP norm (at least in order of magnitude) shift 
from eigenfunctions which localize at a point (zonal spherical harmonics) to eigenfunctions 
which localize on a closed geodesic (highest weight spherical harmonics). It is possible that 
maximal LP growth in the small range of p is related to existence of the same stable elliptic 
geodesies 7 which gave rise to the L p estimates in [TZ2| , p| . 



1. Quantum integrable systems 
We begin by recalling the definition of a completely integrable system. 
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1.1. Completely integrable systems. By a completely integrable system on T*M we 
mean a set of n independent, C°° functions pi, ■ ■ ■ ,p n , on T*M satisfying: 

• {PiiPj} — f° r an 1 < h j < n i 

• dpi A dp2 A • • • A dp n 7^ on an open dense subset of T*M. 
The associated moment map is defined by 

V = ( Pl ,...,p n ):T*M^BcR n . (7) 

We refer to to the set B as the 'image of the moment map.' The Hamiltonians generate an 
action of IR n defined by 

$ t = exp ttE Pl o exp t 2 E P2 • - • o exp t n E Pn . 

We often denote $£-orbits by M. n ■ (x,£). The isotropy group of (x,£) will be denoted by 
1( X £)- When W 1 ■ (x, £) is a compact Lagrangean orbit, then Ti x ^\ is a lattice of full rank in 
]R n , and is known as the 'period lattice', since it consists of the 'times' T G W 1 such that 

$t|a«)(&) = Id. 

We will need the following: 

Definition 1.1. We say that: 

• b G -Bsmg ifP^ib) is a singular level of the moment map, i.e. if there exists a point 
(x,£) G V^ip) with dpi A • • • A dp n (x, £) = 0. S'uc/i a poini (x, £) is called a singular 
point ofV. 

• a connected component ofV~ l (b) (b G B sing )is a singular component if it contains a 
singular point ; 

• an orbit R™ • (x, £) of $ 4 is singular if it is non-Lagrangean, i. e. has dimension < n; 

• b G 5 reg and that P _1 (o) zs a regular level if all points (x,£) G are regular, i.e. 
if dp i A • • • A dp n (x, £) 7^ 0. 

• a component ofV^ 1 ^) ( b G -B s j ns U -B re9y ) regular if it contains no singular points. 

By the Liouville- Arnold theorem [AM], the orbits of the joint flow $ t are diffeomorphic 
to M fc x T m for some (k, m), k + m < n. By the properness assumption on V, a regular level 
has the form 

V-\b) = A (1) (6) U • • • U A M (b), (6 G £ reg ) (8) 

where each A® (6) ~ T™ is an n-dimensional Lagrangian torus. The classical (or geometric) 
multiplicity function m c i(b) = #'P~ 1 (o), i.e. the number of orbits on the level set "P _1 (6), is 
constant on connected components of B reg and the moment map fl7|) is a fibration over each 
component with fiber (§). In sufficiently small neighbourhoods Vt^{b) of each component 
torus, A® (6), the Liouville- Arnold theorem also gives the existence of local action-angle 
variables (if \ I®, 9± , 0^) in terms of which the joint flow of 5 pi , 5 pn is linearized 
[AM]. For convenience, we henceforth normalize the action variables J x , 1$ so that ij = 
0; j = 1, n on the torus A^(6). 

Quantum integrability involves semiclassical pseudo-differential operators, so let us pause 
to set up the notation. Given an open U C IR n , we say that a(x, £; ft) G C°°(?7 x W 1 ) is in 
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the symbol class S m ' k (U x R"), provided 

We say that a G S™' (U x R n ) provided there exists an asymptotic expansion: 

oo 

a(x, £; K) ~ tT m aj(x, 

j=0 

valid for |£| > ~ > with aj(x,C,) G S 0,k ~ j (U x R n ) on this set. We denote the associated 
h Kohn-Nirenberg quantization by Opn(a), where this operator has Schwartz kernel given 
locally by the formula: 

Op h (a)(x,y) = (27th)- n f e^-^' h a(x,& h) d£. 

By using a partition of unity one constructs a corresponding class, Opfi(S m,k ), of properly- 
supported ^-pseudodifferential operators acting on C°°(M). 

Definition 1.2. We say that the operators Qj G Oph(S™ ,k ); j = l,...,n, generate a semi- 
classical quantum completely integrable system if 

[Q t ,Q J ]=0; Vl<i,j<n, 

and the respective semiclassical principal symbols qi,...,q n generate a classical integrable 
system with dq 1 A dq 2 A • • • A dq n ^ on a dense open set fl C T*M — 0. We also assume 
that the finiteness condition (W is satisfied. 

Let us pause to relate this notion to the possibly more familiar notion of n commuting 
(ordinary) pseudo-differential operators Pi, ... , P n . Such operators have a joint spectrum. In 
eigenvalue asymptotics for a joint spectrum, it is natural to consider ladders of eigenvalues: 

{(Ai fc , A nfe ) G Spec(P 1 ,...,P n ); Vj = l,..,n, Km ^ = bj}, (9) 

where |Afe| := J \\ k + ... + A^ fc . We then introduce the semiclassical parameter h with values 
in a sequence {h^ 1 ; k = 1, 2, 3, ...} with %^ — > as k — > oo. Now, consider the semiclassically 
scaled operators: 

Qj- ''I'/- J = l,2,...,n. (10) 

In the semiclassical regime, it is convenient to work with the Q/s rather than the Pj's. For 
the eigenfunctions (f>b,h associated with the ladder @ we clearly have: 

Qj(fib,h = bj(p bth + Hjk4>b,h, where /i jk = o(l) as ft -> 0. (11) 

The operators Qj] j = l,...,n are semiclassical pseudodifferential operators. 

Consequently, in the homogeneous case, the operators P\ = yA, P 2 , P n are classical 
pseudodifferential operators of order one, it readily follows that Qj := HPj G Op(S'°/ 1 ) 
generate a semiclassical quantum integrable system in the sense of Definition |1 . 2| . 
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2. QUASIMODES AND BlRKHOFF NORMAL FORMS 

First, some notation. Following the convention in ||CP|| , we denote ft-microlocal equivalence 
on an open set, Q C T*M by =q. Also, recall that ||UP|| a family of distributions, u^, where 
h G (0, ^o] is called admissible if for any cutoff function, £) G C^°(fl), there is an N G Z 
and sfZ, such that 

\\h N Ph ( X )u h \\ H - s = 0(l) 

uniformly for h G (0, ho]. 

2.1. Inverse Birkhoff normal form. All of our results on localization of modes and quasi- 
modes are based on the use of Birkhoff normal forms. We begin by defining such normal 
forms. 

Let b G B reg and consider the regular level set "P -1 (6) of the moment map with C"' denote 
the pull-back of the Maslov line bundle over the affine torus, T n , given by if = ... = I® = 
and parametrized by the angle variables d\ , ...,9$. We denote the space of smooth sections 
of this bundle by C°°(T n ; £«) . The image of £1® (b) under the map to normal form will be 
denoted by Q =rx5 1) where Di denotes a sufficiently small neighbourhood of G M. n . 
We start with a quantum normal-form construction which can be found in a somewhat 
different form in ||CV2|| . 

LEMMA 2.1. For I = 1, k and j = 1, ...,n, there exist h-Fourier integral operators, : 
C°°(M) — > C°°(T n ; C^), microlocally elliptic on , together with C°° symbols, ff\x\ h) ~ 
EfcLo/ifc(^)ft fc , with f j0 (0) = such that: 

U^*ff\Q 1 -b 1 ,...,Q n -b ni h)U^ =a p ~ 

Moreover, when Qi,...,Q n are formally self- adjoint, the operator C/ 6 can be taken to be 
microlocally unitary. 

Proof. To simplify the writing a little, we drop the superscript (7) in the following and denote 
the h principal symbols of Qi, ...,Q n by q±, ...,q n . Let k : Q(b) — > T n x D\ be a canonical 
map with the property that, 

qjO k(0,I) = qj(I), j = l,...,n, 

for some cjj G C°°(Di). Then, since dq\ , . . . , dq n are linearly-independent on fl(b), by the 
inverse function theorem, there exists fjo G C 00 (R' 1 ) with fja(0) = such that: 

fjo(qi ~ h, ...,q n - b n ) o k = Ij, j = l,...,n. (12) 

By the semiclassical Egorov Theorem, there exists an h-Fouher integral operator Uq : 
L 2 (T n ; C) — > L 2 (M) quantizing k such that for j = 1, ...,n, 

f/(TVi0(Ql-&l,-,Qn-&n)^0=T"xD 1 Ph {Ij ) + H Ph (n (9, I) ) . (13) 

Recall that L 2 (T n ; C) has a natural Hilbert basis given by the shifted exponentials, e m (x) = 
exp(i(m + v/A, x)), with m G Z n . The next step involves reduction of the order of the error 
term on the RHS of (0). We do this by making Fourier series expansions of all symbols in 
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terms of the shifted exponentials, e m (x), and conjugating the expression on the RHS of (fL3l) 
by U x := exp(ihO Ph (v x )) where v x (9,I) G C°°(T n x D x ): 

expl-ihOpniv^] o [0^(1,-) + HO Ph (r x )) o exv[ihOp h {v x )) (14) 
=n [Id - ihOp h {v x )\ o [Opnilj) + fiOp R (n)] o [Id + ihO Ph {v x )\ + 0(ft 2 ). 

Note that, by the pseudodifferential symbolic calculus, 

[Id - ihO Ph ( Vl )] o [Opnilj) + ftOp R (rO] o [Id + ihO Pn (v x )} + 0(ft 2 ) = T n xDl Opfa) + ftOj^-) + 0(£ 2 ) 



Clearly, to get rid of the Oph(e X j) term on the RHS of fll5|) , we would like to solve the 
following transport equations for v x : 

e Xj (9,I) = {v 1 ,I j }(9,I)+r x (9,I). (16) 

Keeping in mind that all operators are acting on L 2 (T n ; C), it follows that these equations 
can be rewritten in the form: 

— Vl (9,I)-^v x (9,I)+r x (9,I) = ;j = l,...,n. (17) 

However, in order to solve (|I7|), we have to subtract off the zeroth shifted Fourier coefficient, 
fji(I), of r x (9,I), and then match Fourier coefficients on both sides of (|i~7|). Consequently, 
if 

fji(qi - h, ■■■An -b n )oK = f jX (I x , ....,!«), 
it follows from di3|)-([l7D above that for j = 1, ...,n, 

Ui'Uv 1 o [f 0j (Q x - b x , Q n - b n ) + hf Xj (Q x - b x , Q n - b n )\ o U U X = TxDl Op n {Ij) + 0(h 2 ) 

(18 

To finish the proof of the lemma, simply apply induction where the inductive step is exactly 
the same as the argument for the first step given above. □ 



2.2. Quasimodes. In this section, we introduce the quasimodes ipb%( x ) := ^& u h 

Suppose that 4>b,n is a joint eigenfunction of Q x , Q n with Qj(pb,h = (bj+Hj(fi))(t>b,h, where 
Hj(h) = o(l). Then, as a consequence of the normal form in Lemma ( |2.1| ), Un =ciW U^* <fib,h 
is an admissible, microlocal solution of the system of equations: 

h~ d 

lW Uh =n o L) m ^ ' Uh ' J ' = ' ' "' ,n ' ' ™ 

where, m^h) G h{Z + § [v®]) + 0(h°°) and m^h) = o(l) as h -> 0. 

Proposition 2.2. The space of admissible microlocal solutions of the system of equations 
in ftTQ ) is one-dimensional. Consequently, the functions 



x) 



u h {9) = exp[im x (K)9 x + ... + im n (h)9 n )} 
span over <C(h) the space of admissible, microlocal solutions to (fJ 
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Proof: Let un, be an admissible, microlocal solution to (|19|) , by using Fourier transforms 
and cutting-off in the frequency variables ||CP|| , one readily shows that there exists a global 
distributional solution, G V ! (T n ;C), to (p!9|) with =q u^. But then, by a simple 
quadrature, we have that vn = c{h) exp[i(mi(h))9i + ... + +m n {h)6 n )} + 0(h°°) for some 
c(h) G C(h). □ 
Proposition ( |2.2| ) has a natural corollary that characterizes the space of the so-called 
semiclassical quasimodes [CV2] attached to a component torus, il"'(b): 

Corollary 2.3. The space of admissible, microlocal solutions to the system of equations 

Qj<f>b,n =n(0 (bj + ^jk{K))<Pb,h + 0(h°°); j = 1, n 
is spanned over C(h) by the functions 

i>l%(x) := uSn(x), 

where, U[ f , is the h-Fourier integral operator constructed in Lemma ( \2. 1\ ) which depends 
locally regularly on b G B reg and un{0) = exp[i(mi(h)9i + ... + m ri {h)O n )}, with rrij{h) G 

z + f [i/(0] + o(h°°). 

Proof: The corollary follows immediately from Proposition ( f2.2| ) together with the Birkhoff 



normal form construction in Lemma ( |2.1| ). □ 
Remark: By applying the lemma of stationary phase in the explicit expressions for the 
functions U^Uf^x) on a sufficiently small open set V G K n , it follows that there exist locally- 
defined phase functions, (pi(x,^) G C°° (V x M. N ), parametrizing the Lagrangian torus, (6) , 
and classical symbols bi(x, £; h) ~ J2j bji(x,^)h 3 , such that: 

4 l) h (x) = (2nh)- N ' 2 f e iMx ^ /h k(x^; h) d£. 



These oscillatory integrals are precisely the usual semiclassical Lagrangian distributions (or 
"quasimodes") attached to Lagrangian tori. 

We still have to determine the asymptotic formulas for the eigenvalues, bj + [ij{K) ■ These 
formulas are usually referred to as the (regular) Bohr-Sommerfeld quantization rules associ- 
ated with a component Lagrangian torus: 

3. Resonant tori and localization of modes 



The main results (Lemmas ( |3.1|) and ( |3.2| ) of this section concern the localization of modes 
on tori. When the fiber 

p-\b) = A (1) (6) U • ■ ■ U A (m ^(6), (20) 

where the A^(b); I = 1, m are has multiple components, it is only true in general that the 
joint eigenfunctions <p\ with approximately the joint eigenvalue A concentrate on the union 
of the tori. The tunnelling problem which concerns us is the following: 

Problem How is the mass of <p\ distributed among the components AS l \b) of V~ l {b)l 

In the case of a symmetric double well potential V(x) = V(—x) on R, it is well-known 
that eigenfunctions are either even or odd, and hence have their mass equally distributed in 
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both wells. A similar phenomenon occurs in any integrable system with symmetries. In this 
case, the wells are resonant, i.e. the germs of the Hamiltonian flow at the two tori (or wells) 
are symplectically equivalent. 

Resonance causes trouble when comparing sup norms of modes and quasi-modes, and 
hence when relating sup norms of modes to the geometry of invariant tori. However, this 
is an exceptional situation: we prove that in generic (non-resonant) situations, modes are 
well-approximated by individual quasimodes and localize on individual tori. To deal with 
the remaining resonant cases, we prove a result (Lemma ( f4.ip ) which allows us to transfer 
attention from modes to quasimodes in the proof of Theorem ( p.l| ). Namely, we prove that 
uniform boundedness of sup norms of modes implies the same for quasimodes ( Lemma 



Much of this material is implicitly known (it appears in a slightly different form in [ |CV2 
|CP|1 in the one-dimensional case), but we have not been able to find the results we need 
explicitly stated and proved. 

3.0.1. Bohr-Sommerfeld, Birkhoff normal form and resonant tori. The notion of resonant 
tori is based on the Bohr-Sommerfeld quantization rules | UV2 | for the joint spectrum of 



Qi, ■■■,Q n - In terms of microlocal analysis, these rules are encoded in the regular quantum 
Birkhoff 'inverse' normal form construction in Lemma ( [2.1|) . Fix a neighbourhood, f2^(6), 
of the torus A^(&) and let (pb,h be a microlocal solution of the system of equations: 

Qj<f>b,h =n(0 ( b j + Vj(ti))<f>b,h+ 0{h°°), with n^h) = o(l) ; j = 1, ...,n. (21) 

Here, by a slight abuse of notation, we have written fij in place of fijk- In addition, suppose 
that 



> CH m , for some m E Z. (22) 



Then, from the 'inverse' Birkhoff normal form construction in Lemma 4.1, it follows that, 
for j = 1, ...,n, 



E /jtW), Vnim*- 1 = 2vrZ + -[!/«] + 0{tf°). (23) 



Let 7^ ; j = 1, n denote the homology generators of the torus T n and k denote the subprin- 
cipal form associated with Qx, Q n (see, for instance [GS2]). Then, to first approximation, 
the equation in (|23|) just the well-known formula: 



The asymptotic identities in (|23|) are called the Bohr-Sommerfeld (B-S) quantization rules for 
the joint eigenvalues (b\ +fii(fi), b n -\-ji n (h)) of (Q\, Q n ) associated with the component 
torus, AW(6 1 ,...,6 n )- 

Note that in Lemma (|2.1|) , we could have just as easily derived the Birkhoff normal form 

as: 

V^giHDe,, HD 6n ; h)V = m Q j - b,, (24) 
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where, gj(xi, x n ; h) ~ J2k 9jk(%i, x n )h? , the g^t's being the inverse functions of the 
fjk's in Lemma fl2.1|) and V : C°°(0^) — >■ x S 1 ) a microlocally unitary ft-Fourier 

integral operator. This way around, we get the (inverse) Bohr-Sommerfeld quantization rules 
expressing the 61 + fa, b n + fjL n in terms of the quantum numbers, rrij(h) G h(Z + n/2[v]) + 
0(h°°): 

00 

EffgW»), ...,m n (^))^ = ^(fi) + 0(/r>) ; j = 1, ...,n. (25) 

In the following, we will need to consider two separate cases, roughly corresponding to 
whether or not the functions fS; I = 1, k (or equivalently, the g^; I — 1, k) are all the 
same for the different component Lagrangian tori. 

We can now define resonance between tori on the same level of the moment map. Let 
b G M n be a fixed regular value of the moment map V and consider the ladder of eigenvalues 

E b (h) := {(61 + /J,i(h), ...A + fa(h)) G Spec(Qx, Q n ); Vj = 1, ...,nfij(h) = o(l)}. 

We say that the component tori, A® (61, 6 n ); Z = 1, are not m resonance if there 
exist eigenvalues (Ai,...,A n ) G Sh(Zz) which satisfy the Bohr-Sommerfeld (B-S) quantization 
rules (p^) associated with the torus, A^(&i, ...,6 n ) to 0(/?. 00 ), but only solve the (B-S) rules 
associated with the other component tori, A^(6i, ...,b n ), with Z 7^ 1 up to some finite order 
m < 00. If a pair of tori on the same level satisfy precisely the same (B-S) quantization 
conditions, then they are said to be resonant. 

3.1. Localization in the non-resonant case. The relationship between quasimodes and 
modes can be fairly complicated when there are resonances between component Lagrangian 
tori. However, as we now show, when the generic non-resonance condition is satisfied, the 
relationship is quite simple. In the following, we let Ob denote the usual order symbol O 
with locally regular dependence on b G B reg . 

Lemma 3.1. Suppose <pb,h is an L 2 -normalized joint eigenfunction of Qi, Q n with joint 
eigenvalues (Ai, A n ) G Hb(h);j = 1, ...,n that satisfies the (B-S) rules to infinite order only 
for the single component torus, A^(6i, b n ) (ie. the component tori, A®(b);l = l,...,m 
are not in resonance). Then, 

(Op h (x {l) )^Op h ( x (l) )hfi) = 1 + o b (h°°). 

Here, Opft(x ) denotes the pseudodifferential cut-off with symbol X ( x i£) supported in 

n^(b). 

Proof: As we have already shown in Proposition ( |2.2|) , the space of microlocal solutions to 
(|19|) is at most one-dimensional. Suppose then that 4>b t n is a joint eigenfunction of Q\, Q n 
with eigenvalues {pi + fa, b n + ji n ) G T,(h). In particular, 4>b,n is then a microlocal solution 
near each component torus; that is, 

Qj(j)b,h =n(Q (h + fa(h))<Pb,h + O b (h°°). 
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Consequently, given Proposition ( |2.2|) , there are two possibilities: In the first case, given 
X®( x :€) e a cutoff supported near the torus, A®, we have that (Optj,(x®)<l>b,hi OPh(x^)4>b,h) > 

Ch m for some m G Z. In such a situation, we already know that the eigenvalues (bi + 
fix, b n + /i n ) G must satisfy the (B-S) quantization rules corresponding to the com- 

ponent torus, A^, to Ob(h°°). The only other possibility is that 

Given the definition of non-resonance above, it follows that the latter scenario holds for I ^ 1 
and the lemma follows. □ 
As a direct consequence of Lemma 3.1 , it follows that in the in non-resonant case, L 2 - 
normalized quasimodes attached to individual Lagrangian tori agree with actual joint-eigenfunctions, 
modulo 0(h°°) errors in L 2 . Moreover, we can now also prove the following correspondence 
principle for the case of lR n actions: 

Proposition 3.2. Let q(x,£) G C£°(T*M) and suppose that the component tori, A^(6); I = 
l,...,m are not in resonance. Then, for all L 2 -normalized joint eigenf unctions, <j>b,n> with 
joint eigenvalues (b\ + fJ>i(h), b n + Unifi)) G T,b(h) with fij(h) = 0{h), there exists a 
component torus, A^(fe), such that: 



(Pp R (q)<f>b,n,(k,h) = W n ! qdfi b + O b (H). 

J\W( b ) 

Here, dfib denotes Liouville measure on h.^ l \b). 



Proof: The proof is similar to the case of torus actions [ZT| . For the sake of completeness, 
we sketch the argument. First, in light of Lemma ( |3.1| ), 

(Opn(q)4>b,n, M = (Op h (q) o Op h ( X {l) )fo,n, Op h { X {l) )h,n) + 0{h°°). (26) 

A standard application of averaging, combined with the semiclassical Egorov theorem and 
the normal form in Lemma gives: 

(Op h (q) X (1) <l>b,H, X {X) 4>b,n) = (27r)" n f qdfi + e(H)+0(h), (27) 

JAW 

where, 

e(h) = (Op h (r)u h ,u h ) 

and r G C °°(T n x Di) with r(6,I) = 0(\I\) and u H {0) = exp[i(mA + ... + m n 9 n )}. An 
integration by parts in the Ii, ...,I n variables shows that: 

(Op h (r)u h ,u h ) = 0{h). 

The proposition follows. □ 



3.2. Localization in the resonant case. In resonant cases, modes do not necessarily 
localize on individual tori. We introduce the language of 'quantum limit measures' to describe 
localization in more general cases. 
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3.2.1. Quantum limits of modes and quasimodes. Let us denote by Mj the set of measures 
on S*M which are invariant under the Hamiltonian M n -action. Among such measures are the 
orbital averaging measures \i Zy defined by J s *m fd^z = liniT-»oo f^r I m ax\tj\<T I* 1 
the case of compact (torus) orbits, one gets the normalized (probability) Lebesgue measure 
on the torus. 

By the 'quantum limit' measures Q of the quantum integrable system, we mean the set 
of weak* limits of the measures d&\ defined by 



(Op(a)(j) X , A ) = / ad$ A , (28) 

JS*M 

where Op(a) is the pseudo differential operator (in some fixed quantization) with complete 
symbol a. It is an easy consequence of Egorov's theorem that Q C Mj. For background, 
terminology and references in a context closely related to this one, we refer to ||TZ|| . 

When b is a regular value of the moment map V, V~ l {b) is given by (|20|), where the 
A^ l '(b);l = l,...,m are n-dimensional Lagrangian tori. Let dfi^u)^ denote the normalize 
Lesbegue measures on such tori. If {(px^ is a sequence of eigenfunctions with p^j} — > b, then 



we have (see WW) 



d§x 3 -> c i^A0')(6), for some c 5 > 0, ^ cj = 1. 

3=1 

It is difficult to determine the coefficients in general. As we have seen, in the generic ('non- 
resonant') case, one can find for each component torus a sequence of eigenfunctions for which 
Cj = 0,j>2. 

Quite analogous to the quantum limit measures are the analogous measures defined for 
quasi-modes associated to regular tori. We thus define the set Q qm C Mj of 'quasimode 
limit measures' to be the weak* limits of the analogous measures d^\ corresponding to 
quasimodes {^a} associated to regular tori. The following lemmas will be crucial for the 
proof of the Theorem in the M. n case: 

Lemma 3.3. For any component Lagrangian torus A® (6), there exists a sequence ip^' n of 
quasimodes such that 

(Op h (q)4% 4%) = (270- L <? d/i b + ChiH), (29) 

JAW 

i.e. that Q qm contains all normalized Lebesgue measures along compact Lagrangean orbits. 

Proof: Since the quasimodes are, by construction, microlocally concentrated on a 
particular component torus, A® (b), Lemma (|3.3j ) is proved in the same way as Proposition 
~T2fo above. □ 

4. Proof of Theorem ( p.l|) 



4.1. Pointwise bounds on quasimodes. In order to prove Theorems p.l| in the possibly 
resonant case, we will need to establish some less direct relations between modes and quasi- 
modes. The next Lemma shows that uniform bounded of sup-norms of modes implies the 
same for quasimodes: 
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Lemma 4.1. Assume ||0a||oo < A for all A e Spec(Pi, P n ). Then for all regular b e R n , 
the L 2 -normalized quasimodes satisfy: 

limsup llrfllU < A m d (b). 

4.1.1. Proof of Lemma Proof: 

Given the component tori A^(b) where / = 1, ...,m, recall that we have denoted by 
the quasimodes in Corollary (|2.3j ) with approximate eigenvalues {b\ + fix, ...,b n + ji n ) e 
which are microlocally concentrated near the component torus, A"' (ft), . By making a 
Gram-Schmidt orthogonalization, we can assume that these functions are orthonormal in 
L 2 . Recall, they satisfy for each j = 1, n: 

in L 2 . Given the (B-S) rules, we know that for any component torus A^(fo), there exist 

e£ > such that for h sufficiently small, the gaps in the (B-S) lattice spectrum are greater 

(i; 

b 



than e^p H. Define 



(0 

e b := mm e\\ 

l<m c i(b) 



and consider the cube, 



C b fi ■= I| i~ e b h + bj + Hj(H), e b h + bj + fij{h)}. 
j=i 

Let Q := (Qi, ...,Q n ) act on L 2 (M) © ■ ■ ■ © L 2 (M) componentwise. Then, we denote by 
Qiffi the (at most) one- dimensional vector space generated by the quasimodes attached to 
A^(6) with quasieigenvalues in Cb,n- In addition, we define 

m c i(b) 

Q b ,n:= Q%- (30) 
l=i 

So, put another way, Q b h is the vector space of functions generated by quasimodes ip^ n ; I = 
1, ...,m c i from all of the component tori with quasieigenvalues in the cube C^n- Similarily, 
we denote by H b ,h the vector space generated by joint eigenfunctions of Q associated with 
joint eigenvalues of Q lying in C^n- As a consequence of the (B-S) rules, we know that 
there are gaps in the joint spectrum of the Q: that is, there exists Cb > such that for all 
fte(0,fio(6)], 

dist ( Spec(Q) — Spec(Q fi C^h) > Spec(Q ) ) > — h. 

Cb 

Let Efc b h denote the spectral projector onto the eigenvalues of Q lying in TZ^n and consider 
a fixed quasimode, ipb,n G Q b ,n- It follows that: 

Q(A,H- En bth A,n) = O b (h°°). (31) 

So, by the spectral theorem, 

ll^ft- E {fo,h,M<kfi\\(o) = °h( Ti0 °)- ( 32 ) 

<t>b,h^Rb,h 



NORMS OF MODES AND QUASI-MODES REVISITED 17 

Conversely, given <pb,h £ 7Zb,h, it follows from the microlocal characterization result in Propo- 
sition that: 



Ubfi- E (<Pb,n,Ar^bM\(0) = O b (h 00 ). (33) 

Let d(Qbfi,lZb,h) = ||n<2 M — ng^I!^,. ||( ) denote the non-symmetric distance between the 
vector spaces Qb,h and lZ b ,n- It follows from (|32|) and (|33|) that, in particular, d(Q bfl , TZ bi n) < 1 
and d(7Zb t h, Qb,h) < 1- Consequently [H], Q bj n — Tib,h- Since Qb t n has dimension uniformly 
bounded independent of h, so does lZb,h- In particular, we have that 

dim TZ b ,h < m d (b). (34) 

Given the estimates in ([H])-(|4]), we have by the Garding and Sobolev inequalities, 

Halloo < II^^Slloo + - - ^ b> ^S)ll(o)- (35) 

The lemma follows from the following two estimates: 

II^SlU = II E^fi^l < A m cl (b) 

(36) 

|| (I - AY (4% - E Rb ^g) II (o) = O b (h°°) 

□ 

4.2. Completion of the proof. The completion of the proof is now exactly as in |[l'Z|| , so 
we only briefly recall it. We first observe: 

Lemma 4.2. Let {ir*dfi A } denote the set of projections to M of normalized Lebesgue mea- 
sures on compact Lagrangean tori A C S*M. Then the family is uniformly bounded as linear 
functionals on L 1 (M). 

Proof: By Lemma ( |3.3|) , each dfi\ is a weak limit of quasi- mode measures corresponding to 
a sequence, say {V'a./i}, of quasimodes. By Lemma ([44]) we have, for V G C(M), 
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oo 



\S M Vit*dn K \ = \im h ^ \(Vip Ajh ,ip A:h ) \ < ||V|| L ilimsup^ H^aaI 

(37) 

<A 2 m c/ (P(A)) 2 \\V\\ L u 

□ 



Corollary 4.3. The invariant Lagrangean tori project without singularities to the base M. 

Proof: Indeed, we have n^dfiA = f\dvol for some fx G L°°(M). But it is straightforward to 
see that if 7t|a has singularities, then fx will blow up along them. □ 

The geometric problem then arises, which completely integrable systems have the property 
that all invariant tori project without singularities to the base Ml It follows immediately 
that M must be a torus, and that no leaf of the Liouville foliation is singular. By a theorem 
of Mane, it then follows that the metric has no conjugate points. But the only metrics on a 
torus without conjugate points are the flat ones (Burago-Ivanov). For more details, we refer 

to Eg. 
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5. Examples of quantum completely integrable systems 

To illustrate some of the variety of quantum completely integrable systems, we describe the 
Liouville foliations of two infinite dimensional families of two-dimensional examples: surfaces 
of revolution and Liouville tori. 

5.1. Spheres of revolution. Suppose that g is a metric on S 2 for which there exists an 
effective action of S 1 by isometries of (S* 2 , g). The two fixed points will be denoted N, S and 
(r, 9) will denote geodesic polar coordinates centered at N, with 9 = some fixed meridian 
7m from N to S. The metric may then be written in the form g = dr 2 + a(r) 2 d9 2 where 
a : [0, L\ — > R + is defined by a(r) = ^\S r (N)\, with |SV(iV")| the length of the distance circle 
of radius r centered at N. For any smooth surface of revolution, the function a satisfies 
a( 2p )(0) = a( 2p \L) = 0,a'(0) = l,a'(L) = —1 and two such surfaces Mi,M 2 are isometric if 
and only if L x = L 2 and ai(r) = a 2 (r) or ai(r) = a 2 (L — r) . 

Complete integrability of H — (i.e. of the geodesic flow) is classical, and follows from 
the existence of the Clairaut integral po(v) : = (v, Since the Poisson bracket {pg, |£| 9 } = 
0, the geodesies are constrained to lie on the level sets of pg; and since both \£\ g and pg are 
homogeneous of degree one, the behaviour of the geodesic flow is determined by its restriction 
to S*S 2 = {\£\ g = 1}. The moment map is given by 

P=m g ,pg):T*S 2 ^R 2 . 

To describe the inverse images, i.e. the invariant tori, we need to specify the function a(r). 
We will assume it is a Morse function on [0, L\ with critical points at {m k } C (0, L). 

5.1.1. Case 1: generic case. For simplicity and with no real loss of generality, let us assume 
that there are two local maxima Mi,M 2 with a(M 2 ) = C 2 > a(Mi) = C\ and one local 
minimum m 2 with a{m 2 ) = c 2 . We note that the canonical involution (x,£) — * (x,—£) 
preserves V~ 1 {p\,p 2 ), and that only the meridian torus P _1 (pi,0) is invariant; otherwise 
there are always at least two components on each level. By homogeneity it suffices to 
consider the level sets P _1 (l, c) in the unit tangent bundle. We will denote the components 
by M c ( ±,fe ) where ± denotes some ordering of each pair of tori interchanged by the involution. 
The following table summarizes the invariant tori, their projections to M and the type of 
singularity of the projection. 



Range of c 

(0) c> C 2 or c < 

(a) c = C 2 

(b) d < c < C 2 

(c) Ci < c < d 

(d) c = ci 

(e) < c < ci 

(f) c = 



components M^'^ 

none 

2 circles 

2 torii 

4 tori 

4 cylinders 

2 tori 

1 (meridian) torus 



projections to M singularities 

Project to closed geodesies singular leaf 

2 horizontal annuli fold 

each projects diffeo to M none 

each projects over all of M none 
projects to 2 horizontal annuli fold 

covers all of M blow-down. 
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5.1.2. Simple surfaces of revolution. Simple surfaces of revolution are examples of toric in- 
tegrable systems [|CV1|| They are the special surfaces where: 

(i) a has precisely one non-degenerate critical point r Q £ (0, L), with a"{r ) < 0, correspond- 
ing to an 'equatorial geodesic' 

(i) the (non-linear) Poincare map V lE for 7^ is of twist type (cf. §1). 



It is proved in [|CV1|1 that the geodesic flow of a simple surface of revolution is is toric 
integrable. 

With the extra assumption on a, the level sets are compact and the only critical level is 
that of the equatorial geodesies 7^ C S*S 2 (traversed with either orientation). The other 
level sets consist of two-dimensional torii. 

As above, let 

(1.1.2) P = (\Z\ g ,pe) : T*S 2 ^B:= {(h, b 2 ) : \b 2 \ < a{r )b x } Clxl + 

be the moment map of the Hamiltonian Reaction defined by the geodesic flow and by 
rotation. The singular set of P is the closed conic set Z : = {(r D , 9, 0,pe) '■ £ [0, 2n),p e £ R}, 
i.e. Z is the cone thru the equatorial geodesic (in either orientation). The image of Z is 
the boundary of B; the map P\r*s g s 2 -z is a trivial S 1 x S 1 bundle over the open convex 
cone B Q (the interior of B). For each h £ B Q , let Hi(Fb,Z) denote the homology of the 
fiber Fb := P _1 (6). This lattice bundle is trivial since B is contractible, so there exists a 
smoothly varying homology basis {71(6), 72(6)} £ Hi(F b ,Z) which equals the unit cocircle 
S^S 2 together with the fixed closed meridian jm when b is on the center line R + ■ (1,0). 
The action variables are given by ||CV1]| , §6, 



1 rr+(b) 



b 2 



(1.1.3) h(b) = f tdx = p e , h(b) = I idx = - f ™ . b\ - -^dr + \b 2 \ 

J nQ>) J l2 {b) 7r Jr-{b) ^ air)' 

where r±(b) are the extremal values of r on the annulus 7r(i*&) (with tt : S*S 2 — > S 2 the 
standard projection). On the torus of meridians in S*S 2 , the value of I 2 equals - and it 
equals one on the equatorial geodesic. So extended, Ji, I 2 are smooth homogeneous functions 
of degree 1 on T*S 2 , and generate 27r-periodic Hamilton flows. It follows that the pair 
X := (Ii,I 2 ) generates a global Hamiltonian torus (S 1 x S' 1 )-action commuting with the 
geodesic flow. We refer to it as the 'standard' Hamiltonian T 2 action on T*S 2 — 0. Up to 
symplectic equivalence, it is unique. 

The singular set of X equals Z := {I 2 = ±pe}, corresponding to the equatorial geodesies. 
The map 

(1.1.4) X : T*S 2 - Z -> T := {(x, y) G R x R + : |x| < y} 

is a trivial torus fibration. Henceforth we often write Tj for the torus X _1 (J) with / £ r o and 
let T = C\Y be the closure of Y as a convex cone. The symplectically dual angle variables 
(<pi = 6, <p 2 ) then give, by definition, the flow times (mod 2n) along the orbits of I 2 ) from 
a fixed point on Fj,, which we may take to be the unique point lying above the intersection 
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of the equator and the fixed meridian on Fb with the geodesic pointing into the northern 
hemisphere. 

The following table summarizes the invariant tori and their projections 

Range of c components M^ ±,k ^ projections to M singularity 

(0) c > C or c < none 

(a) c = C 2 circles Project to closed geodesies singular leaf 

(b) < c < C 2 torii 2 horizontal annuli fold 

(c) c = 1 (meridian) torus covers all of M blow-down. 



5.2. Tori of revolution. Here, we consider a metric g on T 2 which possesses a free effective 
action Rg of S 1 by isometries. We denote the generator of this action by Z. Associated to 
Z is the Clairaut integral pz '■ T*T 2 — > R defined by pziy) = (u, Z). The S l action on T 2 
lifts to T*(T 2 ) as the Hamiltonian flow of pg. Since = 0, we have a Hamiltonian 

G = R 1 x S 1 action ¥ = G" x Rg on T*(T 2 ) with moment map V = {\C\g,Pz)- 
Consider now the level set 

M = V-\l, 0) = {(x, : |£| s = 1, (£, d/dO) = 0} 

It is an affine Lagrangean submanifold of T*(T 2 ) and since it is compact it must be a 
Lagrangean torus. It follows that the stabilizer G m of a point of m G M must be a subgroup 
of the form LZ x {id}. In particular, the G*-orbit of each m G M must be closed. 

As above, let us restrict the natural projection ir : S*T 2 — > T 2 to M, to get tim : M — > M. 
Clearly, 7r M is equivariant with respect to the S 1 action Rg on T 2 and its lift. It sends the 
G*-orbits of m G M to closed geodesies of T 2 and since pz = on M, these geodesies are 
orthogonal to the -R^-orbits, which we will call 'parallels'. Further, we note that the tangent 
space T m M is spanned by the Hamilton vector fields , resp. E pz . Since itm* sends these 
vectors to orthogonal vectors, itm is a submersion and hence a covering map. Its degree is 
one because of the equivariance and the fact that it sends G*-orbits to closed geodesies (this 
also follows from Lalonde's theorem cited above). Hence it is a diffeomorphism. We will 
refer to the geodesies in M as the meridians. 

Fixing a point m G M and x Q = 7Tm("t. ) G T 2 we get a coordinate system on T 2 with 
coordinates (r,9) : x = (r, 6) <^=> <&( r ' e \m ) = ir]^~(x). Since the meridians are orthogonal 
to the 'parallels', this forms an orthogonal coordinate system for g. Hence g has the form: 
g = dr 2 + a(r) 2 d9 2 . 



5.2.1. Liouville torii ||B.K.S||[|KMSfl . A Riemannian torus (M,g) of dimension 2 is called a 
Liouville torus if the metric has the form 

g = [Ui(xi) - U 2 (x2)](dxl + dx\) 

in the coordinates (xi,x 2 ) on [0, 1] x [0, 1]. The functions Uj are assumed to be periodic of 
period 1 and to satisfy Ui(x 2 ) — ^2(^2) > 0. The geodesic flow is the Hamilton flow on T*M 
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generated by 



(^ + £ 2 2 ). 



[Ui(xx)-U 2 (x 2 ) 



e 2 . 



[U 1 (x 1 )-U 2 (x 2 ) 
It Poisson commutes with the first integral 

We then have the homogeneous moment map 

V = (pi >P2 ) = {H, ^):T*M-0^ R 2 . 

We note that the canonical involution (x,£) — > preserves V~ 1 (pi,p2), and that no 

torus is invariant; hence there are always at least two components on each level. As above, 
it suffices to consider the level sets in the unit tangent bundle, which we continue to denote 

by M^' k l 

To describe the inverse images, i.e. the invariant tori, we need to specify the functions 
Uj(xj). We will assume that Uj is a Morse function on [0, 1] with critical points at {cjk} C 
(0, 1) with critical values {vjk}- We briefly review the description of the tori under the 

and one minimum Cj_ and that 
We thus have 



assumption that each Uj has precisely one maximum Cj 

[|KMS 



the critical points do not coincide (see 
Ci+ > Cn > c 2+ > c 2 _. 

Range of c 

(0) c > ci + or c < c 2 _ 

(a) c = ci+ 

(b) ci_ < c < c 1+ 

(c) c = Ci_ 

(d) c 2+ < c < Ci_ 
(a') c = c 2+ 
(b') c 2 - < c < c 2+ 
(c') c = c 2 _ 



components M^'^ 
none 
2 circles 
2 torii 
4 cylinders 
4 tori 

analogous to (a) 
analogous to (b) 
analogous to (c) 



for a detailed description), 
projections to M 

Project to closed geodesies 
2 horizontal annuli 
each projects over all of M 
each projects diffeo to M 
2 circles 

2 vertical annuli 



singularity 

singular leaf 
folds 
none 
none 

singular leaves 
folds 



5.2.2. Tori of revolution: Quantum theory. Consider the metric on M 2 /Z 2 

g = a(x)(dx 2 + <i£ 2 ) 

where a is a periodic function. Translations in the £ variable are then isometries so this is a 
torus of revolution. The Laplacian 

A = -A 

a 

commutes with <9/<9£ and hence there is an ONB of joint eigenf unctions of the form 

4 NtX = e 2 * iN ^ N}X (x). 
The eigenvalue problem is separable in these variables and gives 

ip" - Att 2 N 2 ^) = -Xaip. 
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Bourgain [Bour] observes that for a = 1 — Tsin 2 7rx and for sufficiently large N there are 
eigenfunctions satisfying 

||0a||oo>A5, ||0 A || i2 = l. 

These eigenfunctions concentrate near the projection of a stable geodesic given by the 
parallel x — 0. Note that Theorem 0.2 gives in this case a lower bound of ||0a||oo > C(e)\*~ e 
for any e > 0. Since da(x) = — 27rsin7rx cos7ra; and since x G [0,1] there are two critical 
circles: at x = and at x — |, a maximal parallel and a minimal parallel (which are 
geodesies). There are four embedded cylinders consisting of geodesies which spiral in towards 
the latter circle. 



5.2.3. Liouville tori: Quantum theory [B.K.S|1[|KMS|1 . The Laplace Beltrami operator on a 
Liouville torus, Q, with metric g = [Ui(xi) — U 2 (x 2 )](dxi + dx^) is just: 



A=[U 1 (x 1 )-U 2 (x 2 )]- 1 (—,+ 



8T 
dxo 



The commuting operator with principal symbol S(x,£) in 5.4 is just: 



S = [U 1 (x 1 ) - U 2 {x 2 )Y l U 2 (x 2 )-^ + U ± (xi 



2 



dx{ dx 

The joint eigenfunctions of A and S can be described by two coupled O.D.E.(see K.M.S.). 
Moreover, in the case where c = C\ + or c = c 2+ the two circles in the joint level set are 
stable and elliptic and project to a single geodesic on Q. In the case c = Ci + this projected 
geodesic is identified with the line given by Xi = Mi where U\(xi) attains its non-degenerate 
maximum at this point, whereas, when c = c 2+ this geodesic gets identified with the line 
x 2 = M 2 where U 2 (x 2 ) attains its maximum at M 2 . There are eigenfunctions attached to 
these levels which satisfy ||0a||oo > As and again this is consistent with the prediction given 
by Theorem 0.2. Note also that there exist hyperbolic geodesies on the level sets given by 
c = Ci_ and c = c 2 - which are limit sets for the geodesies on these cylinders. These are the 
lines given by X\ = mi and x 2 = m 2 respectively, where mi and m 2 are the points where 
Ui and U 2 achieve their respective minima. In this case, [Tl], one can show that there are 
joint eigenfunctions concentrating near the projections of these hyperbolic orbits and satisfy 
|| 0a || oo > As(logA)~2 ; slightly better than the ||0a||oo > C(e)As~ e given by Theorem 0.2. 



5.3. Other examples. Other classical examples include the Euler top and geodesic on 
asymmetric ellipsoids in any dimension. These examples are also quantum completely inte- 
grable. Tri-axial ellipsoids have special umbilic points where all directions are geodesic loop 
directions. Hence they are candidates for maximal eigenfunction growth. The loops form 
a Lagrangean manifold with boundary which has a blow-down projection over the umbilic 
points. We do not believe however that the eigenfunctions actually do have maximal sup 
norm growth. 



NORMS OF MODES AND QUASI-MODES REVISITED 



23 



6. Further discussion and open problems 



In [[TZjl , we list a number of (still) open problems in this area. We would like to mention 



three of them here since they seem rather glaringly open. 

• Are the results on quantum completely integrable systems valid under the weaker as- 
sumption of classical complete integrability? E.g. if A is a quantum integrable Lapla- 
cian, and if its eigenfunctions have a given growth rate, can the additition of a lower 
order term such as A + V with V G C°°(M) change the growth rate? 

• Is the existence of a stable elliptic closed geodesic sufficient for the existence of a se- 

n— 1 

quence {0a} of eigenfunctions satisfying | |0a| |oo ~ A - ^ - ? Note that if would be sufficient 
to obtain sup norm estimates for KAM systems. 

• How do the results generalize to Schroedinger operators A + V on R n ? Here, we may 
assume V grows at infinity to ensure that the spectrum is discrete. 
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